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Embedded  Ultrasonic  Structural  Radar  with  Piezoelectric  Wafer 
Active  Sensors  for  Damage  Detection  in  Cylindrical  Shell  Structures 
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Structural  health  monitoring  (SHM)  is  a  major  component  of  the  vehicle  health  monitoring  (VHM)  concept 
currently  considered  for  the  civilian  and  military  aerospace  applications.  Piezoelectric  wafer  active  sensors 
(PWAS)  are  one  of  the  candidate  embedded  sensors  considered  for  SHM  applications.  PWAS  are  inexpensive, 
non-intrusive,  un-obtrusive,  devices  that  can  be  used  in  both  active  and  passive  modes.  In  active  mode,  PWAS 
generated  Lamb  waves  that  can  be  used  for  damage  detection  through  pulse-echo  or  pitch-catch  techniques. 

An  efficient  application  of  the  pulse-echo  method  with  PWAS  technology  is  through  the  phased  array 
technique.  In  the  embedded  ultrasonics  structural  radar  (EUSR)  concept,  an  array  of  closely  spaced  PWAS  is 
used  to  detect  structural  cracks  based  on  the  scanning  beam  principle  with  guided  Lamb  waves.  A  scanning 
beam  of  ultrasonic  Lamb  waves  works  like  ^‘structural  radar”.  When  encountering  a  crack,  it  generates 
echoes  and  backscatter.  This  concept  was  initially  developed  by  the  authors  for  guided  Lamb  waves  traveling 
in  flat  plates.  In  this  paper,  the  EUSR  concept  is  extended  to  cylindrical  shells.  First,  the  theory  of  guided 
waves  in  cylindrical  shells,  with  the  associated  modes  and  dispersive  frequencies  is  reviewed.  It  is  shown  that 
cylindrical  shells  accept  three  types  of  guided  ultrasonic  waves:  longitudinal,  flexural,  and  torsional.  The  first 
and  second  of  these  ways  can  be  associated  with  the  Lamb  waves  in  fiat  plates,  while  the  third  can  be  related 
to  the  shear-horizontal  guided  waves  in  fiat  plates.  Subsequently,  the  paper  describes  validation  experiments 
performed  on  cylindrical  shells  of  various  curvatures.  It  is  shown  that  the  EUSR  concept  works  on  cylindrical 
shells  with  curvatures  representative  to  actual  aircraft  structure  just  as  well  as  it  works  on  fiat  plates. 


I.  Introduction 

Embedded  nondestructive  evaluation  (NDE)  is  an  emerging  technology  that  will  allow  the  transition  of  the 
conventional  ultrasonics  methods  to  embedded  structural  health  monitoring  (SHM)  systems  such  as  those  envisioned 
for  vehicle  health  management  (VHM).  SHM  for  VHM  requires  the  development  of  small,  lightweight,  inexpensive, 
unobtrusive,  minimally  invasive  sensors  to  be  embedded  in  the  airframe  with  minimum  weight  penalty  and  at 
affordable  costs.  Such  sensors  should  be  able  to  scan  the  structure  and  identify  the  presence  of  defects  and  incipient 
damage.  Current  ultrasonic  inspection  of  thin  wall  structures  (e.g.,  aircraft  shells,  storage  tanks,  large  pipes,  etc.)  is  a 
time  consuming  operation  that  requires  meticulous  through-the-thickness  C-scans  over  large  areas.  One  method  to 
increase  the  efficiency  of  thin-wall  structures  inspection  is  to  utilize  guided  waves  (e.g..  Lamb  waves)  instead  of  the 
conventional  pressure  waves.  Guided  waves  propagate  along  the  mid- surface  of  thin-wall  plates  and  shallow  shells. 
They  can  travel  at  relatively  large  distances  with  very  little  amplitude  loss  and  offer  the  advantage  of  large-area 
coverage  with  a  minimum  of  installed  sensors.  Guided  Lamb  waves  have  opened  new  opportunities  for  cost-effective 
detection  of  damage  in  aircraft  structures,  and  a  large  number  of  papers  have  recently  been  published  on  this  subject. 
Traditionally,  guided  waves  have  been  generated  by  impinging  the  plate  obliquely  with  a  tone-burst  from  a  relatively 
large  ultrasonic  transducer.  Snell’s  law  ensures  mode  conversion  at  the  interface,  hence  a  combination  of  pressure 
and  shear  waves  are  simultaneously  generated  into  the  thin  plate.  However,  conventional  Lamb-wave  probes  (wedge 
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and  comb  transducers)  are  relatively  too  heavy  and  expensive  to  consider  for  widespread  deployment  on  an  aircraft 
structure  as  part  of  a  SHM  system.  Hence,  a  different  type  of  sensors  than  the  conventional  ultrasonic  transducers  is 
required  for  the  SHM  systems.  Piezoelectric  wafer  active  sensors  (PWAS)  are  inexpensive,  non-intrusive,  un¬ 
obtrusive  devices  that  can  be  used  in  both  active  and  passive  modes.  In  active  mode,  PWAS  generated  Lamb  waves 
that  can  be  used  for  damage  detection  through  pulse-echo  or  pitch-catch  techniques. 


Figure  1  Reference  coordinates  and  characteristic  dimensions  of  a  hollow  cylinder  in  which  guided  waves  are  propagated 
(after  ref.  1 0) 


II.  Guided  Waves  in  Thin-wall  Cylindrical  Shells 

In  flat  plates,  ultrasonic  guided  waves  travel  as  Lamb  waves  and  as  shear  horizontal  (SH)  waves.  Lamb  waves  are 
vertically  polarized,  while  SH  waves  are  horizontally  polarized.  Lamb  waves  can  be  symmetrical  or  antisymmetrical 
with  respect  to  the  plate  midplane.  At  lower  frequency-thickness  products,  only  two  Lamb  wave  types  exist:  So, 
which  is  a  symmetrical  Lamb  wave  type  resembling  the  longitudinal  waves;  and  Aq,  which  is  an  antisymmetric  Lamb 
wave  type  resembling  the  flexural  waves.  At  higher  frequency-thickness  values,  a  number  of  Lamb  waves  are 
present,  S„  and  A„,  where  ^  =  0,  1,  2,  ...  At  very  high  frequencies,  the  So  and  Ao  Lamb  waves  coalesce  into  the 
Rayleigh  waves,  which  are  confined  to  the  plate  upper  and  lower  surfaces.  The  Lamb  waves  are  highly  dispersive 
(wave  speed  varies  with  frequency);  however.  So  waves  at  low  frequency-thickness  values  show  very  little 
dispersion.  The  SH  waves  are  also  dispersive,  with  the  exception  of  the  first  mode,  SHo,  which  does  not  show  any 
dispersion  at  all.  Ultrasonic  guided  waves  in  flat  plates  were  first  described  by  Lamb^  A  comprehensive  analysis  of 
Lamb  wave  was  given  by  Viktorov^,  Achenbach^,  Graff^,  Rose^,  and  Royer  &  Dieulesaint^. 

The  guided  Lamb  waves  have  the  important  property  that  they  stay  confined  inside  the  walls  of  a  thin-wall 
structure,  and  hence  can  travel  over  large  distances.  In  addition,  the  guided  waves  can  also  travel  inside  curved  walls; 
this  property  makes  them  ideal  for  applications  in  the  ultrasonic  inspection  of  aircraft,  missiles,  pressure  vessel,  oil 
tanks,  pipelines,  etc.  The  study  of  guided  waves  propagation  in  cylindrical  shells  can  be  considered  a  limiting  case  of 
the  study  of  guided  waves  propagation  in  hollow  cylinders.  As  the  wall  thickness  of  the  hollow  cylinder  decrease 
with  respect  to  its  radius,  the  hollow  cylinder  approaches  the  case  of  a  thin-wall  cylindrical  shell.  Several 
investigators  have  considered  the  propagation  of  waves  in  solid  and  hollow  cylinders.  Love^  studied  wave 
propagation  in  an  isotropic  solid  cylinder  and  showed  that  three  types  of  solutions  are  possible:  (1)  longitudinal;  (2) 
flexural;  and  (3)  torsional.  At  high  frequencies,  each  of  these  solutions  is  multimodal  and  dispersive.  Meitzler^ 
showed  that,  under  certain  conditions,  mode  coupling  can  exist  between  the  various  wave  types  propagating  in  solid 
cylinders  such  as  wires.  Extensive  numerical  simulation  and  experimental  testing  of  these  phenomena  was  done  by 
Zemenek^. 

Comprehensive  work  on  wave  propagation  in  hollow  circular  cylinders  was  done  by  Gazis^^.  A  comprehensive 
analytical  investigation  was  complemented  by  numerical  studies.  The  nonlinear  algebraic  equations  and  the 
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corresponding  numerical  solutions  of  the  wave-speed  dispersion  curves  were  obtained.  These  results  work  found 
important  applications  in  the  ultrasonic  nondestructive  evaluation  (NDE)  of  tubing  and  pipes.  Silk  and  Bainton^^ 
found  equivalences  between  the  ultrasonic  in  hollow  cylinders  and  the  Lamb  waves  in  flat  plates  and  used  them  to 
detect  cracks  in  heat  exchanger  tubing.  Rose  et  al.^^  used  guided  pipe  waves  to  find  cracks  in  nuclear  steam  generator 
tubing.  Alleyne  et.  al.^^  used  guided  waves  to  detect  cracks  and  corrosion  in  chemical  plant  pipework. 

A  brief  review  of  the  mathematical  modeling  and  the  main  results  for  guided  waves  in  cylindrical  shells  are  given 
next^^’^^’  The  coordinates  and  characteristic  dimensions  are  shown  in  Figure  1;  a  and  b  are  the  inner  and  outer  radii 
of  the  tube,  and  h  is  the  tube  thickness.  The  variables  r,  6,  and  z  are  the  radial,  circumferential,  and  longitudinal 
coordinates.  The  modeling  starts  from  the  equation  of  motion  for  an  isotropic  elastic  medium,  in  invariant  form: 


//Au  +  (yl-}-//)VV«u  =  pid^uldt^) 


(1) 


where  u  is  the  displacement  vector,  p  is  the  density,  X  and  p  are  Lame’s  constants,  and  A  is  the  three-dimensional 
Laplace  operator.  The  vector  u  is  expressed  in  terms  of  the  dilation  scalar  potential  (j)  and  the  equivolume  vector 
potential  H  according  to: 

u  =  V^  +  VxH  (2) 

V.H  =  F(r,0  (3) 

For  free  motion,  the  displacement  equations  of  motion  are  satisfied  if  the  potentials  ^  and  H  satisfy  the  wave 
equations 

(4) 

(5) 


cIaU- 


where  Cp  =(/l  +  2//)/ p  and  c\=  p! p  are  the  pressure  and  shear  wave  speeds,  respectively.  We  express  the 
potentials  and  the  wave  equations  in  cylindrical  coordinates  as 

(j)-  f{r)Qo^n0QO^{Q)t  +  ^z) 

(r)  sin  nO  +  ^z) 

Hq  =  ge  {r)  cos  nO  +  ^z) 

=  g^  (r)  sin  nO  co^{cot  +  ^z) 

where  wave  motion  along  the  z  axis  with  wave  number  ^  is  assumed.  Substitution  of  Equation  (6)  into  Equations  (4) 
and  (5)  yields 


(6) 


{ls.  +  Qr ! Cp)(j)  =  Q 
(A  +  ry^/cs)//,  =0 


(7) 


+0}^  lcl)H^-{2lr^){dHglde)  =  0 
(A-l/r^  +CO^  lcl)Hg+{2lr^){dH^  Idd)  =  0 
where  A  =  is  the  Laplace  operator.  The  following  notations  are  introduced: 

=0?  I  (8) 

The  general  solution  is  expressed  in  terms  of  Bessel  functions  J  and  Y  or  the  modified  Bessel  functions  /  and  K  of 
arguments  a^r  =  \ar\  and  =  \Pr\  as  determined  by  Equation  (8)  are  real  or  imaginary.  The  general  solution  is  of 
the  form 

g,=A,-ZMr)  +  B,WMr) 
g,=Hgr-go)  =  A-  {P,r)  +  B,  ■  (Ar) 
g2=i{gr+ge)  =  A- 2„_i (Pp)  +  B,  ■  iP,r) 

where,  for  brevity,  Z  denotes  a  J  or  /  Bessel  function,  and  W  denotes  sl  Y  or  slK  Bessel  function,  as  appropriate.  Thus, 
the  potentials  are  expressed  in  terms  of  the  unknowns  A,  B,  A],  Bj,  A2,  B2,  A3,  B3.  Two  of  these  unknowns  are 
eliminating  using  the  gauge  invariance  property  of  the  equivolume  potentials.  The  strain-displacement  and  stress- 
strain  relations  are  used  to  express  the  stresses  in  terms  of  the  potential  functions.  Then,  we  impose  the  free-motion 
boundary  conditions 

(7^^  =  (7^^  =  (7^^  =  0  at  r^a  and  r  =  b  (10) 


(9) 
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Thus,  one  arrives  at  a  linear  system  of  six  homogeneous  equations  in  six  unknowns.  For  nontrivial  solution,  the 
system  determinant  must  vanish,  i.e., 

|c..|  =  0,/j=l  ...6  (11) 

The  coefficients  Cy  in  Equation  (11)  have  complicated  algebraic  expressions  (Gazis,  1959)  that  are  not  reproduced 
here  for  sake  of  brevity.  We  retain,  however,  that  a  characteristic  equation  exists  in  the  form: 

Q.^{a,b,X,jLi,fd,c)  =  ^  (12) 

where  a  and  b  represent  the  inner  and  outer  radii  of  the  tube,  while  X  and  ju  represent  the  Lame  constants.  This 
implicit  transcendental  equation  is  solved  numerically  to  determine  the  permissible  guided-wave  solutions.  As  shown 
by  Gazis^^,  three  basic  families  of  guided  waves  exist: 

•  Longitudinal  axially  symmetric  modes,  L(0,  m),  m  =  1,  2,  3,  . . . 

•  Torsional  axially  symmetric  modes,  T(0,  m) ,  m  =  1,  2,  3,  . . . 

•  Flexural  non-axially  symmetric  modes,  F(^,  m),  n  =1,2,  3,...  m  =  l,2,  3,... 

Within  each  family,  an  infinite  number  of  modes  exist  such  that  their  phase  velocities,  c,  for  a  given  frequency- 
thickness  product,  fd,  represent  permissible  solutions  of  an  implicit  transcendental  Equation  (12).  The  index  m 
represents  the  number  of  the  mode  shape  across  the  wall  of  the  tube.  The  index  n  determines  the  manner  in  which  the 
fields  generated  by  the  guided  wave  modes  vary  with  angular  coordinate  0  in  the  cross-section  of  the  cylinder.  For 
the  F  modes  of  family  n,  each  field  component  can  be  considered  to  vary  as  either  sin(^^  or  cos(^^.  It  is  also 
observed  that  the  index  n  represents  the  mode  shape  of  flexing  of  the  tube  as  a  whole. 

The  longitudinal  and  torsional  modes  are  also  referred  to  as  the  axial  symmetric  or  ^  =  0  modes.  These  axial 
symmetric  modes  are  preferred  for  defect  detection  in  long  pipes  since  the  pipe  circumference  is  uniformly 
insonified.  The  longitudinal  modes  are  easier  to  generate  with  conventional  ultrasonic  transducers.  They  are  good  for 
the  location  of  circumferential  cracks.  However,  for  the  location  of  axial  cracks  and  corrosion,  the  torsional  modes, 
though  more  difficult  to  generate  with  conventional  ultrasonic  transducers,  are  recommended. 

Thus,  the  existence  in  thin-wall  cylinders  of  guided  waves  similar  to  the  Lamb  waves  present  in  flat  plates  was 
identified.  An  examination  of  the  differential  equations  dependence  on  the  ratios  h/r  and  h/X  indicates  that,  for 
shallow  shells  Qi/r  and  h/X  «  1),  the  longitudinal  modes  approach  the  Lamb  wave  modes,  whereas  the  torsional 
modes  approach  the  SH  modes.  In  fact,  it  can  be  shown  that  the  L(0,I)  mode  corresponds  to  the  Aq  Lamb  mode, 
while  the  L(0,2)  mode  corresponds  to  the  So  Lamb  mode,  etc.  The  flexural  modes  remain  a  mode  type  specific  to 
tube  waves  and  without  equivalence  in  flat  plate  waves. 


Figure  2:  (a)  conventional  Lamb  wave  ultrasonic  transducer  (~25  grams,  ~$400,  40-mm  tall);  (b)  piezoelectric  wafer  active 
sensors  (PWAS):  0.08  grams,  $15,  0.2  mm  thin. 
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III.  Guided-Wave  Generation/Detection  with  Piezoelectric  Wafer  Active  Sensors 


In  conventional  ultrasonics,  guided  waves  are  generated  by  impinging  the  structural  surface  obliquely  with  a 
ultrasonic  beam  from  a  relatively  large  ultrasonic  transducer  affixed  to  a  wedge.  Snell’s  law  ensures  mode  conversion 
at  the  interface,  hence  a  combination  of  pressure  and  shear  waves  are  simultaneously  generated  into  the  structure.  If 
the  structure  is  thin  walled,  guided  plate  of  shell  waves  are  being  created.  Another  method  of  creating  guided  waves 
in  a  structure  is  by  using  a  comb  transducer.  The  comb  spacing  is  such  that  it  tunes  with  the  guided  wave  half 
wavelength.  However,  conventional  Lamb-wave  probes  (wedge  and  comb  transducers)  are  relatively  too  heavy  and 
expensive  to  be  considered  for  widespread  deployment  on  an  aircraft  structure  as  part  of  a  SHM  system.  A  typical 
wedge  transducer  weighs  25  grams,  costs  ~$400,  and  is  40-mm  tall  (Figure  2a).  Hence,  a  different  type  of  sensors 
than  the  conventional  ultrasonic  transducers  is  required  for  the  SHM  systems. 

In  recent  years,  piezoelectric  wafers  permanently  attached  to  the  structure  have  been  extensively  used  for  the 
generation  and  detection  of  guided  waves.  These  simple  devices  are  inexpensive,  very  lightweight,  and  unobtrusive 
(Figure  2b).  They  seem  ideally  suited  for  structural  health  monitoring  applications.  They  can  act  as  both  transmitters 
and  receivers  of  ultrasonic  waves.  They  are  called  piezoelectric  wafer  active  sensors  (PWAS).  PWAS  weight  around 
68  mg,  are  0.2  mm  thick,  and  cost  only  ~$15  each. 
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(a) 

Figure  3  PWAS  interaction  with  So  and  Ao  Lamb  modes 


(b) 


PWAS  operated  on  the  piezoelectric  principle  that  couples  the  electrical  and  mechanical  variables  in  the  material 
(mechanical  strain,  mechanical  stress,  Tu,  electrical  field,  and  electrical  displacement  Df  in  the  form: 


D, 


^kij^k 
~  ^jkl’^kl  ^jk^k 


(13) 


E  f 

where  is  the  mechanical  compliance  of  the  material  measured  at  zero  electric  field  (E  =  0),  is  the  dielectric 

permittivity  measured  at  zero  mechanical  stress  (T  =  0),  and  dky  represents  the  piezoelectric  coupling  effect.  For 
embedded  NDE  applications,  PWAS  couple  their  in-plane  motion,  excited  by  the  applied  oscillatory  voltage  through 
the  piezoelectric  effect,  with  the  Lamb-waves  particle  motion  on  the  material  surface.  Lamb  waves  can  be  either 
quasi-axial  (SO,  SI,  S2,  ...),  or  quasi-flexural  (AO,  Al,  A2,  ...)  as  shown  in  Figure  3.  PWAS  transducers  can  act  as 
both  exciters  and  detectors  of  the  elastic  Lamb  waves  traveling  in  the  material.  PWAS  can  be  used  as  both  active  and 
passive  probes  and  thus  can  address  four  SHM  needs: 

1 .  Active  sensing  of  far-field  damage  using  pulse-echo,  pitch-catch,  and  phased-array  methods 

2.  Active  sensing  of  near-field  damage  using  high-frequency  impedance  method 

3.  Passive  sensing  of  crack  initiation  and  growth  through  acoustic  emission 

4.  Passive  sensing  of  damage-generating  events  through  detection  of  low-velocity  impacts 

An  important  characteristic  of  PWAS,  which  distinguishes  them  from  conventional  ultrasonic  transducers,  is  their 
capability  of  tuning  into  various  guided  wave  modes.  Though  a  complete  modeling  of  this  interaction  is  not  yet 
available,  some  simplified  models  exist  that  clarify  the  underlying  principles.  Giurgiutiu^"^  developed  a  1-D  analysis 
of  the  PWAS-structure  interaction  using  the  space-domain  Fourier  analysis.  This  model  illustrates  the  principles  of 
PWAS  Lamb-wave  mode  tuning,  and  opens  the  path  for  a  more  comprehensive  analysis.  The  analysis  starts  with 
transformed  equations: 
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dy^  dy^  cl  cl 

where  ^  and  y/  are  the  transformed  potential  functions,  cl  =  (A-\-2ju)/ p  and  cl-pIp  are  the  longitudinal 
(pressure)  and  transverse  (shear)  wave  speeds,  X  and  p  are  Lame  constants,  p  is  the  mass  density.  The  excitation  is 
assumed  in  the  form  of  time-harmonic  shear  stress,  applied  by  the  PWAS  to  the  structural  surface  over  the 

interval  -a  <  x  <  +a.  Assuming  perfect  bonding  between  PWAS  and  structure  concentrates  the  shear  transfer  effects 
at  the  PWAS  ends.  Hence,  the  shear  stress  ^x)  and  its  Fourier  transform  f(^)  are  given  by 

t(x)  =  uTq  [S(x-a)-S(x-\-a)] ,  f(^)  =  uTq  [-2/sin^a]  (15) 

where  ^(x)  is  the  Dirac  function.  Equations  (14),  (15)  accept  the  Fourier-domain  solution: 
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where  is  the  axial  strain,  is  the  axial  displacement,  and 


) cos ph cos qh,  D^={^  -q  y  cos phsmqh  +  A^  pqsm ph cos qh 


(16) 


(17) 


Na  -  '^h)  sin  ph  sin  qh,  D^=(^  -q  y  sin  ph  cos  qh  +  4^  pq  cos  ph  sin  qh 

Applying  the  inverse  Fourier  transform  in  the  space  domain,  and  adding  the  time-domain  harmonic  behavior  yields 
the  forward  wave  solution: 
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The  integral  in  Equation  (18)  is  singular  at  the  roots  of  Ds  and  Da.  These  roots  correspond  to  the  symmetric  and  anti¬ 
symmetric  Lamb  wave  modes  having  eigenvalues,  ,^2  ,^2  The  integral  in  Equation  (18)  is 

evaluated  through  the  residue  theorem,  using  a  semicircular  contour  consisting  in  the  upper  half  of  the  complex  ^ 
plane.  Hence, 
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Figured  (a)  Predicted  Lamb  wave  strain  amplitude  in  a  1-mm  aluminum  plate  under  a  7-mm  PWAS  excitation;  (b) 
experimental  verification  of  excitation  sweet  spot  at  300  kHz  (ref.  14) 


The  summations  in  Equations  (19)  cover  all  the  symmetric  )  and  anti-symmetric  (^^)  Lamb  wave  modes  that 
exist  for  a  given  value  of  in  a  given  plate^"^.  A  plot  of  this  solution  in  the  0-1000  kHz  bandwidth  is  presented  in 
Figure  4a.  The  sin^a  contained  in  Equation  (19)  displays  maxima  when  the  PWAS  length  4  =  2a  equals  an  odd 
multiple  of  the  half  wavelength,  and  minima  when  it  equals  an  even  multiple  of  the  half  wavelength.  Several 
such  maxima  and  minima  exist,  each  associated  with  a  certain  Lamb  mode  of  its  own  wavelength.  These  minima  and 
maxima  allows  us  to  achieve  Lamb-mode  tuning.  We  can  tune  into  a  selected  Lamb  mode  by  choosing  the 
appropriate  frequency  and  PWAS  geometry.  Figure  4a  illustrates  the  tuning  of  the  So  mode  at  300  kHz.  This  happens 
because  at  this  frequency,  the  amplitude  of  the  Aq  mode  goes  through  zero  while  that  of  the  So  is  still  strong.  Thus, 
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we  have  tuning  of  the  So  mode,  and  rejection  of  the  Aq  mode.  Figure  4b  shows  the  experimental  confirmation  of 
this  So  “sweet  spot”.  Lamb  wave  mode  tuning  offers  considerable  advantages.  It  allows  us  to  select  Lamb  wave 
modes  that  are  most  appropriate  for  the  particular  application  being  considered.  The  tuning  of  the  So  mode  is 
important  for  the  detection  of  certain  defects,  such  as  cracks.  The  So  mode  is  less  dispersive  than  the  Ao  mode  at  low 
fd  values  and  produces  better  wave  packets.  When  the  Ao  mode  is  rejected,  the  remaining  So  mode  gives  loud  and 
clear  pulse-echo  reflections  from  the  through-the-thickness  cracks^ 


IV.  Phase  Array  Principles  Used  in  the  Embedded  Ultrasonic  Structural  Radar 


An  important  application  of  the  PWAS  based  ultrasonics  is  that  based  on  phased  array  principles.  Giurgiutiu  and 
Bao^^  developed  a  PWAS  phased  array  application  and  named  it  embedded  ultrasonics  structural  radar  (EUSR).  Its 
principle  of  operation  is  derived  from  two  general 
concepts: 

(1)  The  of  guided  Lamb  wave  generation  with 
piezoelectric  wafer  active  sensors  (PWAS) 

(2)  The  principles  of  conventional  phased-array 
radar. 

The  phased  array  beam  forming  principles  can  be 
easily  understood  through  the  following  3 -element 
array  example  presented  in  Figure  5.  In  this  example, 
three  wave  sources  were  placed  half  wave  length 
apart.  When  activated,  the  waves  generated  by  each 
source  would  interfere  in  the  wave  field,  generating 
constructive  and  destructive  interference  patterns.  The 
beam  in  the  middle  corresponds  to  where  the 
constructive  interference  is  at  a  maximum.  By  timing 
the  relative  phase  between  the  wave  signals  emanating 
from  each  source,  the  beam  of  maximum  constructive 
interference  can  be  steered  around  the  horizon. 

Lamb  waves  can  exist  in  a  number  of  dispersive 
modes.  However,  through  smoothed  tone-burst 
excitation  and  frequency  tuning,  it  is  possible  to 
confine  the  excitation  to  a  particular  Lamb  wave 
mode,  of  carrier  frequency  fc,  wave  speed  c,  and  wave 

length  X-cj  .  Hence,  the  smoothed  tone-burst  signal  generated  by  one  PWAS  is  of  the  form: 


Figure  5 


10  15  20  25  30  35 

Simulated  beam  forming  from  a  3-element  array 


Sj  (^)  =  Sq  (t)  cos  Itt  fj 


(20) 


where  sq  (t)  is  a  short-duration  smoothing  window  that  is  applied  to  the  carrier  signal  of  frequency  fc  between  0  and 


tp. 

The  principles  of  conventional  phased-array  radar^^  are  applied  to  the  PWAS  array,  assuming  a  uniform  linear 
array  of  M  PWAS,  with  each  PWAS  acting  as  a  point-wise  omni-directional  transmitter  and  receiver.  The  PWAS  in 
the  array  are  spaced  at  the  distance  d,  which  is  assumed  much  smaller  than  the  distance  r  to  a  generic,  far-distance 
point,  P.  Since  d  «  r,  the  rays  joining  the  sensors  with  the  point  P  can  be  assimilated  with  a  parallel  fascicle,  of  (j) . 
Therefore,  for  the  PWAS,  the  distance  will  be  shorted  by  m{d  cos^) .  If  all  the  PWAS  are  fired  simultaneously, 
the  signal  from  the  PWAS  will  arrive  at  P  quicker  by  A^(^)  =  m(d  cos^)/c .  Yet,  if  the  PWAS  are  not  fired 
simultaneously,  but  with  some  individual  delays,  ,  m  =  0,1,..., M-1,  then  the  total  signal  received  at  point  P  will 
be: 

1  f  r  ^ 

=  —  +  (21) 

Vr  ^=0  V  J 

where  1/r  represents  the  decrease  in  the  wave  amplitude  due  to  the  omni-directional  2-D  radiation,  and  rjc  is  the 
delay  due  to  the  travel  distance  between  the  reference  PWAS  (m  =  0)  and  the  point  P.  (Here  wave-energy 
conservation,  i.e.,  no  dissipation,  is  assumed.) 
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(22) 


Transmitter  beamforming:  if  we  have  =  mA(^) ,  then  Equation  (21)  becomes: 

s^{t)  =M  I  ^““1 

Vr  'v  c) 

That’s  to  say,  there  is  an  M  times  increase  in  the  signal  strength  with  respect  to  a  simple  sensor.  This  leads  directly  to 
the  beamforming  principle  that  if  =  md cos(^q)/c  ,  and  since  =  md cos(^)/c ,  then  constructive  interference 

(beamforming)  takes  place  when  cos(^)  =  cos(^q)  ,  i.e.  at  angles  ^  ^  •  Thus,  the  forming  of  a  beam  at 

angles  (j)^  and  -(j)^  is  achieved  through  delays  in  the  firing  of  the  sensors  in  the  array. 

Receiver  beamforming:  if  the  point  P  is  an  omni-directional  source  at  azimuth  ,  then  the  signals  received  at  the 
m^^  sensor  will  arrive  quicker  by  mA^^  =  md cos(^o)/c  .  Hence,  we  can  synchronize  the  signals  received  at  all  the 
sensors  by  delaying  them  by: 

=  cos((Z>o)  (23) 

c 
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.  A 1 

Transmission  puise 
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Figure  6  The  basis  of  pulse-echo  method:  (a)  transmitted  smooth-windowed  tone-burst  of  duration  tp]  (b)  received  signal 
to  be  analyzed  for  the  duration  k,  starting  at  tp,  in  order  to  identify  the  time  of  flight  delay,  t. 


Pulse-echo  method:  Assume  that  a  target  exists  at  azimuth  (/>q  and  distance  R.  The  transmitter  beamformer  is 


sweeping  the  range  in  increasing  angles  (j)  and  receives  an  echo  when  The  echo  will  be  received  on  all 

sensors,  but  the  signals  will  not  be  synchronized.  To  synchronize  the  sensors  signals,  the  delays 
Sm  {(l>Q )  =  md  cos(^Q )  /  c  need  to  be  applied. 

The  original  signal  is: 


5p(0 


(24) 


At  the  target,  the  signal  is  backscattered  with  a  backscatter  coefficient,  A.  Hence,  the  signal  received  at  each  sensor 
will  be: 


AM 

R 


2^  A 


The  receiver  beamformer  assembles  the  signals  from  all  the  sensors  with  the  appropriate  time  delays: 

2R 


^  m=0 


t- 


-  +  ^  (d)-S 


(25) 


(26) 


Constructive  interference  between  the  received  signals  is  achieved  when  =  md  cos(^q  )  /  c  .  Thus,  the  assembled 
receive  signal  will  be  again  boosted  M  times,  with  respect  to  the  individual  sensors: 
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(27) 


(0  ~ 


AM 
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2  M-1  A 

“  X  ‘^T  U 


Then  the  target  range  can  be  calculated  as: 


R  =  ct/2  (28) 

Therefore,  what  we  need  to  determine  from  the  signals  is  the  time  delay  between  the  receive  signal,  ^r(^)  and  the 


transmit  signal  sjQ).  Figure  6  shows  the  transmitted  and  received  signals. 
9-element  PWAS  array 


HP-331 20A 
signal 
Generator 


TDS-210  digital 
oscilloscope 


Data 

acquisition 
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Figure  7  Experimental  setup  for  the  EUSR  experiment 


V.  Experimental  Results 

To  verify  the  theoretical  results,  experiments  were  conducted  at  the  University  of  South  Carolina  in  the  Laboratory 
for  Adaptive  Materials  and  Smart  Structures  (LAMSS).  The  experimental  setup  is  presented  in  Figure  7.  The 
experimental  setup  has  been  identical  with  that  used  previously  in  the  development  of  the  EUSR  concept  on  a  flat 
plate.  It  consists  of  a  1-mm  thick  aluminum  plate  of  square  geometry  (1220  mm  by  1200  mm).  A  PWAS  array  was 
constructed  in  the  middle  of  the  plate.  The  PWAS  array  consisted  of  nine  piezoelectric  wafers  (APC  850 
piezoceramic,  0.2  mm  thick,  7  mm  by  7  mm  square).  The  PWAS  array  was  wired  to  a  connection  pad.  An  HP- 
33 120A  signal  generator  was  used  to  produced  a  3-count  smoothed  tone  burst  of  300  kHz.  The  excitation  frequency 
was  such  as  to  tune  into  the  SO  Lamb  wave  mode.  A  TDS-210  digital  oscilloscope  was  used  to  collect  and  digitize  the 
signals.  The  signal  generation  and  collection  proceeded  in  a  round-robin  fashion.  In  turns,  one  PWAS  was  transmitter 
and  all  the  other  PWAS  were  receiver.  The  signal  was  collected  on  all  the  receivers  and  also  on  the  transmitter.  Then, 
the  next  PWAS  in  the  array  became  transmitter  and  the  process  repeated.  Thus,  a  total  of  9x9  =  81  elemental  signals 
were  collected.  These  signals  were  processed  with  the  EUSR  algorithm.  Figure  8a  shows  the  location  and  size  of  the 
simulated  crack  (narrow  slit).  Figure  8b  shows  the  imaging  of  the  crack  on  the  EUSR  GUI  and  the  reconstructed 
pulse-echo  signal  (A  scan).  The  signal  reconstruction  and  the  crack  imaging  were  done  in  virtual  time  by  applying 
the  EUSR  algorithm  to  the  set  of  81  elemental  signals.  Figure  8b  shows  the  A  scan  of  the  crack  echo  recorded  as  the 
virtual  beam  points  to  an  azimuth  angle  of  90^.  Loud  and  clear  reception  of  the  echo  signal  is  observed.  During 
interactive  tests,  it  was  found  that  this  echo  diminishes  rapidly  as  the  beam  direction  shifts  to  the  right  and  to  the  left 
of  the  simulated  crack  target.  By  stepping  the  scanning  beam  through  the  interval  0^  to  180^,  a  complete  scan  of  the 
plate  upper  half  could  be  obtained.  Putting  together  all  these  beams  in  a  polar  plot,  and  using  a  color  scale  of  the 
signal  amplitude  we  were  able  to  create  a  mapped  image  of  the  plate  upper  half  In  this  mapped  image,  the  crack 
location  is  clearly  defined  by  a  pronounced  color  change  (Figure  8b). 
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Figure  8  Crack  detection  at  zero  curvature  (flat  plate):  (a)  location  and  size  of  simulated  crack;  (b)  crack  imaging  and 
crack  echo  in  the  EUSR  algorithm 


Table  1  Curvature  of  the  specimen 


R 

c 

AL=L-c 

100.00E+3 

1219.99 

50.00E+3 

1219.97 

20.00E+3 

1219.81 

10.00E+3 

1219.24 

8.80E+3 

1219.02 

1 

5.00E+3 

1216.98 

3.90E+3 

1215.03 

5 

2.75E+3 

1210.02 

10 

2.24E+3 

1204.98 

15 

2.00E+3 

1201.17 

1.95E+3 

1200.20 

20 

1.00E+3 

1145.73 

L  =  1220,  all  in  mm 


In  order  to  verify  the  effect  of  curvature,  the  thin  aluminum  plate  was  temporary  bent  into  cylindrical  shapes  using 
thin  wire  ropes  and  tightening  screws.  The  realized  curvatures  were  related  to  the  shortening  of  the  chord  with 
respect  to  the  original  flat  length.  After  curving,  the  original  flat  length,  L,  became  the  circular  arc,  while  the  shorted 
wire  rope  became  the  chord,  c.  The  difference  between  the  original  length,  L,  and  the  shortened  length,  c,  was 
denoted  by  AL.  The  chord  c  can  be  calculated  with  the  formula 

c  =  27?sin—  (29) 

2R 

Then,  the  radius  of  curvature,  R,  could  be  calculated  with  the  formula: 


R  = 


L-AL 


(30) 
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The  values  of  curvature  obtained  with  Equation  (30)  are  given  in  Table  1.  In  our  experiment  we  went  through  the 
shortenings  AL  =  5mm,  10mm,  15mm  and  20mm  corresponding  to  curvature  values  R  =  3.90  m,  2.75  m,  2.24  m,  1.95 
m.  These  curvatures  were  applied  in  two  directions.  One  direction  was  with  the  chord  parallel  to  the  PWAS  array 
(direction  1),  while  the  other  was  with  the  chord  perpendicular  to  the  PWAS  array  (direction  2). 


Figure  9  Crack  detection  at  high  curvature  (R  =  1 .95  m)  in  direction  1 :  (a)  photo  of  specimen  while  being  bend;  (b)  crack 
imaging  and  crack  echo  in  the  EUSR  algorithm 


VI.  Summary  and  Conclusions 

This  paper  has  discussed  the  use  of  the  embedded  ultrasonics  structural  radar  (EUSR)  with  piezoelectric  wafer 
active  sensors  (PWAS)  for  damage  detection  in  cylindrical  shell  structures.  This  work  is  in  continuation  of  previous 
work  of  similar  nature  performed  on  flat  plates.  In  Section  I,  a  brief  introduction  highlighting  the  advantages  of  using 
guided  waves  for  nondestructive  evaluation  (NDE)  of  thin-wall  structures  was  presented.  Also  discussed  was  the 
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need  for  unobtrusive  low-cost  and  light  weight  transducers  to  be  used  in  embedded  structural  health  monitoring 
(SHM)  system.  The  PWAS  transducers  were  identified  as  strong  candidates  for  performing  this  function. 

Section  II  discussed  guided  waves  in  thin- wall  cylindrical  shells.  The  mathematical  model  for  the  analysis  of 
wave  propagation  in  hollow  cylinders  was  reviewed  and  its  general  solution  was  presented.  It  was  shown  that  the 
general  solution  accepts  three  wave  types:  (i)  longitudinal  axially  symmetric  modes,  L(0,  m),  m  =  1,  2,  3,  ...;  (ii) 
torsional  axially  symmetric  modes,  T(0,  m)  ,  m  =  1,  2,  3,  ...;  and  (iii)  flexural  non-axiahy  symmetric  modes,  ¥(n, 
m),  n  =1,2,  3,  ...,  m  =  l,2,  3,  ...  An  examination  of  the  differential  equations  dependence  on  the  ratios  h/r  and 
h/A  indicates  that,  for  shallow  shells  (h/r  and  h/A  «  1),  the  longitudinal  modes  approach  the  Lamb  wave  modes, 
whereas  the  torsional  modes  approach  the  SH  modes.  In  fact,  it  can  be  shown  that  the  L(0,I)  mode  corresponds  to  the 
Ao  Lamb  mode,  while  the  L(0,2)  mode  corresponds  to  the  So  Lamb  mode,  etc.  The  flexural  modes  remain  a  mode 


Figure  10  Crack  detection  at  high  curvature  (R  =  1.95  m)  in  direction  2:  (a)  photo  of  specimen  while  being  bend;  (b)  crack 
imaging  and  crack  echo  in  the  EUSR  algorithm 


Section  III  discussed  the  guided- wave  generation  and  detection  with  PWAS.  It  was  shown  that  due  to  their 
bidirectional  piezoelectric  mechano-electrical  coupling,  PWAS  can  act  as  both  generators  and  detectors  of  elastic 
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waves.  PWAS  can  perform  four  SHM  functions:  (i)  active  sensing  of  far-field  damage  using  pulse-echo,  pitch-catch, 
and  phased-array  methods;  (ii)  active  sensing  of  near-field  damage  using  high-frequency  impedance  method;  (iii) 
passive  sensing  of  crack  initiation  and  growth  through  acoustic  emission;  (iv)  passive  sensing  of  damage-generating 
events  through  detection  of  low- velocity  impacts.  The  mathematical  model  of  a  PWAS  mounted  to  the  surface  of  a 
flat  plate  was  presented.  The  general  solution  derived  through  the  space-domain  Fourier  transform  was  reviewed. 
For  the  case  of  perfect  bonding  between  the  PWAS  and  the  structure,  a  closed  form  solution  was  obtained.  The 
closed- form  solution  is  expressed  as  a  summation  of  the  individual  responses  of  several  Lamb  wave  modes.  In  this 
summation,  the  individual  contribution  of  each  Lamb  wave  mode  varies  with  frequency.  A  plot  of  this  solution  vs. 
frequency  indicated  that  Lamb  wave  tuning  opportunities  exist  whereas,  for  a  given  PWAS  geometry,  certain 
frequencies  exist  that  maximize  the  response  of  a  particular  Lamb  wave  mode  with  respect  to  the  other  modes.  For  a 
1-mm  thick  aluminum  plate,  it  was  predicted  that  the  high-dispersion  AO  mode  will  be  maximized  at  around  100  kHz 
while  the  low-dispersion  SO  mode  will  be  maximized  at  around  300  kHz.  These  theoretical  predictions  were  verified 
very  well  by  experimental  results. 

Section  IV  treated  the  phased  array  principles  applied  to  arrays  of  PWAS.  The  discussion  was  started  with  the 
analysis  of  the  wave  field  generated  around  three  PWAS  that  were  excited  with  a  slight  phase  difference.  It  was 
shown  that  a  ultrasonic  beam  can  be  formed  through  the  constructive  interference  of  the  wave  fields  generated  by  the 
PWAS.  The  beam  angle  could  be  adjusted  by  changing  the  relative  phase  of  the  three  elementary  signals.  These 
principles  were  extended  to  an  array  of  M  identical  equally-spaced  PWAS.  Using  the  parallel-ray  assumption,  the 
transmission  beamforming  and  reception  beamforming  algorithms  were  presented.  Eventually,  the  formula  for  pulse- 
echo  detection  of  defects  using  a  scanning  beam  principle  was  derived.  These  principles  were  used  to  create  the 
embedded  ultrasonics  structural  radar  (EUSR)  algorithm  with  a  graphical  user  interface  (GUI)  for  virtual-time 
processing  of  phased-array  signals. 

In  Section  V,  experimental  results  were  presented.  A  9-element  PWAS  array  was  constructed  on  a  1-mm  thick 
2024-T3  aluminum  plate.  A  small  narrow  slit  was  manufactured  in  the  plate  to  simulate  crack  damage.  Using  general 
purposes  laboratory  instrumentation,  a  set  of  8 1  elemental  signals  was  collected  through  round-robin  excitation  of  the 
individual  PWAS  and  collection  of  signals  on  all  the  PWAS.  The  signals  corresponded  to  tuned  excitation  of  the  low- 
dispersion  SO  Lamb  wave  mode.  The  elemental  signals  were  processed  in  virtual  time  with  the  EUSR-GUI.  Good 
and  accurate  imaging  of  the  crack  was  obtained.  Then,  the  effect  of  curvature  on  the  crack  detection  capability  was 
studied.  Four  curvature  values  were  tested.  The  radius  of  curvature  for  these  curvatures  was  3.9  m  (13  ft);  2.75  m  (9.2 
ft);  2.24  m  (7.5  ft);  1.95  m  (6.5  ft). The  curvatures  were  applied  in  two  directions:  (1)  with  the  chord  parallel  to  the 
PWAS  array;  (2)  with  the  chord  perpendicular  to  the  PWAS  array.  Good  crack  detection  was  obtained  for  all 
curvatures  and  all  directions.  No  significant  difference  due  to  curvature  could  be  observed. 

We  conclude  that  the  EUSR  method  with  PWAS  phased  arrays  can  be  successfully  applied  to  detect  and  image 
simulated  cracks  in  cylindrical  shell  structures.  Further  work  should  focus  on  refining  the  theoretical  and 
experimental  results  and  then  testing  the  method  on  realistic  aerospace  structures. 
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